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ABSTRACT 


A  mass  accelerator  system  using  a  series  of  imploded  annular  plasma 
pinches  to  propel  projectiles  of  mass  ranging  from  grams  to  kilograms  up 
to  high  velocities  is  analyzed  theoretically.  Such  a  device  would  have 
applications  ranging  from  basic  materials  and  impact  studies  to  space  probe 
launching  and  impact  fusion. 


I.  INTRODUCTION 

The  implosion  dynamics  of  self  z-pinched  annular  discharges  is  a  subject 
that  has  been  well  studied  in  the  context  of  such  devices  as  the  plasma 
focus  ,  or  dynamically  imploded  annular  foils  or  gas  puffs  for  radiation 
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sources.  In  this  paper  we  give  some  analysis  of  a  recently  proposed  concept 
for  a  mass  accelerator  system  using  a  series  of  imploded  discharges  (MAID), 
which  are  axially  aligned  and  imploded  sequentially  on  to  the  surface  of  a 
suitably  tapered  projectile  in  order  to  propel  it.  The  system  appears  capable 
of  accelerating  projectiles  of  mass  ranging  from  grams  to  kilograms  up  to 
velocities  well  above  10®  cm/sec  and  possibly  as  high  as  10^  cm/sec.  It  is 
also  expected  to  have  good  projectile  survival  and  stability  properties, 
could  have  a  high  repetition  rate,  and  the  attainable  accelerations  appear 
to  be  limited  principally  by  the  requirement  that  the  projectile  should  not 
undergo  material  damage  such  as  crushing  or  spalling.  The  potential  applica¬ 
tions  which  would  motivate  the  development  of  this  technology  include  basic 
materials  and  impact  studies,  the  launching  of  small  space  probes,  or  possibly 
impact  fusion. 

The  basic  accelerator  module  is  shown  in  Figure  1.  It  consists  of  an 
anode  and  cathode  which  for  example  could  be  discs  with  a  hole  at  the  center 
of  sufficient  radius  to  allow  passage  of  the  projectile.  At  a  suitable  time 
interval  before  arrival  of  the  projectile  in  the  anode-cathode  gap  from  the 
left  in  Figure  1,  a  high  voltage  pulse  is  switched  on  across  the  electrodes 
which  is  sufficient  to  initiate  a  discharge  through  flashover  in  the  low 
density  gas  along  the  inner  surface  of  the  insulator  which  separates  the 
electrodes.  This  high  current  annular  discharge  is  then  driven  by  its  azi¬ 
muthal  magnetic  field  radially  inward  away  from  the  insulator  surface  and 
accelerates  towards  the  axis  of  the  module.  As  the  discharge  advances 


radially  inward,  it  accumulates  and  compresses  more  plasma  by  sweeping  up 
the  background  gas  via  the  snow  plow  effect,  so  that  the  velocity  and  mass 
of  gas  colliding  with  the  projectile  surface  depends  on  the  initial  back¬ 
ground  gas  pressure  (for  example  <^10  torr)  and  module  size.  The  inner  edge 
of  the  insulator  between  the  anode  and  cathode  could  also  be  slightly  tapered 
as  shown,  so  that  the  imploding  annular  plasma  becomes  itself  more  tapered 
during  its  implosion  so  as  to  match  the  projectile  surface,  although  this 
tapering  is  not  necessary  and  the  insulator  could  also  be  a  simple  straight 
cylindrical  section  separating  the  anode  and  cathode.  The  discharge  is  timed 
so  that  it  collides  with  the  projectile  surface  when  its  front  edge  is  approxi¬ 
mately  aligned  with  the  cathode.  The  projectile  is  propelled  both  by  the  im¬ 
ploded  discharge  momentum  and  pressure,  and  by  reaction  to  its  ablated  surface 
material. 

Figure  2  shows  a  mass  accelerator  consisting  of  a  series  of  axially 
aligned  discharge  modules.  The  voltages  of  each  module  would  be  typically 
kilovolts  to  ten  kilovolts  depending  on  the  application.  As  the  projectile 
is  propelled  forward  from  one  anode-cathode  gap  into  the  next  module,  a 
sensing  device  (for  example  the  interruption  of  a  light  beam)  is  used  to 
determine  when  to  trigger  the  next  diode  modules.  This  is  timed  so  that  the 
next  discharge  arrives  at  the  projectile  surface  when  it  is  appropriately 
positioned  in  the  anode-cathode  gap  for  further  acceleration,  and  the  sequence 
is  repeated  as  the  projectile  advances  axially  along  the  accelerator.  Note 
that  for  devices  that  take  full  advantage  of  the  power  compression  that 
results  if  the  radius  to  the  inner  edge  of  the  insulator  is  much  larger  than 
the  projectile  radius,  several  discharges  ahead  of  the  projectile  could  be 
initi ated  and  be  in  different  stages  of  implosion  at  any  one  time. 
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Some  advantages  of  this  system  are  (i)  the  implosions  effectively  focus 
the  discharge  gas  momentum  onto  the  projectile  where  the  maximum  pressures 
in  the  system  then  result,  (ii)  power  compression  (discussed  in  Section  2), 
(iii)  the  expected  good  survivability  properties  of  either  metal ic  or 
insulator  projectiles  propelled  to  high  velocities  by  gas  pressure  instead 
of  by  eddy-current  inducing  magnetic  fields,  and  (iv)  frictional  contact 
with  confining  walls  as  occurs  in  rail  guns  is  avoided,  so  that  the  system 
has  a  high  potential  repetition  rate. 

As  an  example  of  the  approximate  scaling  involved,  assume  that  the  pro¬ 
jectile  has  a  mass  of  M  grams  and  is  accelerated  up  to  a  velocity  of  U  cm/sec 
with  an  average  acceleration  of  U  cm/sec^  (averaged  over  a  module  plus  gap 
between  modules).  The  average  total  power  in  the  electrical  energy  pulse 
that  is  required  to  generate  this  acceleration  is  then 


P  =  =  10"^  MU  U/e  watts. 


where  e  is  the  overall  efficiency  with  which  electrical  energy  is  transferred 
from  the  electrical  storage  system  to  projectile  kinetic  energy.  In  this 
formula  we  distinguish  between  P  and  the  power  of  each  module,  where 
n  is  the  average  number  of  discharge  modules  that  are  operating  at  any  one 
time  with  their  plasmas  in  various  stages  of  implosion  timed  to  impact  the 
projectile  surface  as  it  passes  through.  This  power  would  produce  a  projec¬ 
tile  with  final  kinetic  energy. 

nA 

S-IO*^  M  (— 1  Joules 
VlOV 


if  applied  for  an  acceleration  length  L  =  U|/20  for  the  case  in  which  the 
acceleration  0  is  constant.  For  example,  assuming  U  -  5. 10®  cm/sec^. 
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£  *  0.1,  =  3.10®  cm/sec  gives  for  the  final  maximum  power  =  1.5.10® 

M/n  watts  per  module,  and  a  projectile  energy  of  4.5.10®  M  Joules  acquired 
in  an  accelerator  length  of  100  meters. 

More  generally  the  efficiency  in  (1)  (and  acceleration)  are  not  constant 
along  the  accelerator.  Further,  e  is  made  up  of  two  parts, 

^  “  ^impl  ^prop  ’ 

where  £.  ,  is  the  efficiency  with  which  electrical  energy  is  transferred  to 

impl 

the  imploding  plasma  shells,  and  Spj,Qp  is  the  propulsion  efficiency  with 
which  this  energy  is  in  turn  transmitted  to  projectile  kinetic  energy.  These 
efficiencies  are  functions  of  various  parameters  characterizing  the  modules 
and  the  projectile,  and  will  be  discussed  in  Sections  II  and  III. 

The  projectiles  (Figure  3)  could  be  constructed  out  of  various  materials. 
For  cases  in  which  high  currents  are  used,  it  is  important  to  minimize  heating 
inside  the  projectile  (if  it  is  electrically  conducting)  which  can  arise  from 
eddy  currents  in  its  outer  layers  due  to  penetration  of  the  surface  magnetic 
fields  through  the  discharge  plasma  into  the  solid  projectile,  since  this 
can  cause  premature  melting  of  the  projectile  before  a  high  velocity  is 
reached.  This  problem  is  alleviated  for  metalic  projectiles  by  arranging 
that  the  plasma  pressure  on  the  projectile  derive  mainly  from  the  implosion 
energy  of  the  discharge  so  that  the  magnetic  fields  remain  small  at  the 
inner  plasma  edge  during  impact  with  the  projectile  surface  (the  driving 
current  could  in  fact  be  zero  at  impact  as  discussed  in  Section 
2D,  which  also  appears  desirable  for  projectile  stability.  Section  IV). 

Eddy  current  heating  can  also  be  eliminated  for  cases  in  which  lower 
accelerations  (weaker  projectiles)  are  acceptable  by  constructing  the 
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projectile  of  insulator  material .  In  this  situation,  current  is  then  able 
to  flow  only  through  the  surfjice  plasma  layer  pinched  against  the  projectile, 
so  that  Joule  dissipation  is  zero  inside  the  projectile  except  for  the  thin- 
conducting  skin  layer  that  has  been  converted  to  plasma  at  its  outer  edge. 

This  plasma  layer  advances  into  the  projectile  relatively  slowly  via  thermal 
conduction  and  ablation  during  the  acceleration,  i.e.,  more  slowly  than  would 
be  the  case  for  the  more  deeply  penetrating  skin  current  layer  that  occurs 
for  metal ic  projectiles  in  devices  in  which  magnetic  fields  provide  the  primary 
propulsion  instead  of  plasma  momentum  and  reaction  forces  as  in  this  device. 

II.  IMPLOSION  RELATIONS 

The  implosion  speed  of  the  annular  discharge  as  it  approaches  the 

projectile  surface  must  satisfy  >  U  tan  0,  where  U  is  the  projectile 

velocity  and  9  its  cone  half-angle.  Under  these  circumstances  the  timing 
can  be  arranged  so  that  contact  is  not  made  between  the  discharge  and  any 
part  of  the  tapered  projectile  prior  to  arrival  of  the  projectile  in  the  gap, 
and  further,  the  discharge  can  be  timed  to  collide  with  at  least  some  of  the 
projectile  surface  before  the  projectile  has  had  time  to  move  ahead  of  the 
electrode  gap.  We  also  note  that  the  implosion  time. 


^impl 


'^impl^®^ 


for  the  plasma  to  implode  radially  from  the  insulator  surface  to  the  projec¬ 
tile  surface,  can  be  larger  than  the  transit  time  of  approximately  A/U  taken 
for  the  projectile  to  cross  the  electrode  gap,  where  the  inwardly  directed 
radial  velocity  of  the  annular  discharge,  V|j^p.|  =  | aj ,  is  a  function  of  its 
radial  position  a,  a  is  the  radius  of  the  insulator  surface,  a  the  radius 

W  r 
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of  the  projectile  at  a  point  located  at  a  distance  ~il/2  back  from  the  leading 
edge  of  the  tapered  section  of  the  projectile,  and  l  the  axial  length  of  the 
discharge  as  shown  in  Figure  1.  For  cases  in  which  this  condition  is  satis¬ 
fied  some  power  multiplication  can  be  obtained  in  that  the  driving  voltage 
and  Current  can  be  applied  so  that  the  discharge  plasma  accumulates  implosion 
momentum  for  a  time  longer  than  the  encounter  time  available  to  communicate 
this  accumulated  momentum  to  the  projectile.  It  is  also  important  that  the 
gas  density  in  which  the  discharge  is  initiated  be  sufficiently  low  that 
the  final  implosion  velocity  exceeds  the  projectile  velocity,  while  at  the 
same  time  sufficiently  high  for  the  implosion  radius  involved  that  the  dis¬ 
charge  momentum  provides  a  high  surface  pressure  on  impact  with  the  projec¬ 
tile.  Analysis  of  these  requirements  involves  solving  hydrodynamic  and 
circuit  equations  for  magnetically  imploded  plasma  shells.  We  next  give  a 
simple  generalization  of  a  model  by  Mosher  for  the  case  of  interest  to  us, 
namely  that  in  which  snow  plowing  of  ambient  gas  occurs. 


A.  Approxi mate  Imp! osi on  Model . 

The  equations  of  motion  for  a  cylindrically  symmetric  gas  implosion 
in  a  module  of  radius  a^  and  gap  length  I  can  be  written 


where 


is  the  integrated  radial  momentum,  the  fluid  velocity. 


(4) 


(5) 


p  =  nT  (1  +  Z)  (6) 
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the  gas  pressure. 


is  the  total  plasma  current. 

If  we  assume  a  snow  plow  model  in  which  all  the  gas  of  initial  uniform 
density  is  swept  up  by  the  implosion,  then  P^Trp^A  (a^  -  a^)a.  Further, 
assuming  the  annular  shell  thickness  Aa  is  small  compared  with  a  (this  breaks 
down  near  the  end  of  the  implosion  and  is  discussed  in  Section  2F),  the 
coupled  system  of  imploding  plasma  and  circuit  equations  become 


(9) 


V  =  I  =  ^  (LI)  +  IR.  I  =  -  Q  (10) 

L  =  L  A)  .  (11) 

g  ^2  a 

R*®i!,/2TTaa6  ,  (12) 
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for  an  implosion  driven  by  discharge  of  a  capacitor  C,  where  L,  Lg,  R, 

V  are  the  total  inductance,  inductance  external  to  the  module,  discharge 
resistance,  and  voltage,  and 

m„  =  Traf,  ilp„  ( 

g  w  0 


is  the  mass  of  gas  in  the  module.  In  the  formula  for  the  discharge  resis¬ 
tance  we  take  6  as  the  current  penetration  skin  depth. 


^2+  % 

=  (^-^) 
'4ira' 


+  6 


flash 


(14) 


where  a  is  assumed  constant  and  is  the  small  initial  conducting  layer 

thickness  associated  with  the  initiating  surface  breakdown  at  a^. 


B .  Efficiency  of  Implosion. 

Multiplying  (10)  by  I  and  integrating  gives 

T  X 

y*dt  VI  =  J  dt  ih  L^I^)  +y*dt  I^R  . 

0  0  0 

Using  I  =  -  Q,  Q  =  CV,  L  =-2ila/ac2  ,  and  B  =  2I/ac  with  l(t=0)  =  0,  we  find 


2C(V^  -  V^) >  %  LI' 


T 

■/ 


dt  (|fp)  2TTaa5, 


T 

./ 


dt  I'R 


(15) 


which  displays  the  partition  of  energy  between  the  change  in  capacitor 
storage,  inductive  energy,  kinetic  energy  of  the  implosion  and  Joule 

dissipation, 
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The  efficiency  for  transfer  of  electrical  energy  to  implosion  kinetic 
energy  will  be  defined  as 


£-• 


impl 


1 


W. 


impl 


(>sCVj)  (1-f) 


^■mpl 


t,. 


impl 

dt  I^R  +  %  ^^impl^- 


(16) 


In  this  expression,  f  is  the  fraction  of  the  original  capacitor  bank  energy 
that  remains  stored  in  the  capacitor  at  the  end  of  the  implosion,  assuming 
a  switch  can  be  opened  at  that  time  so  that  the  energy  can  be  usefully  pre¬ 
served  for  later  use.  For  example,  the  voltage  may  reverse  prior  to  impact 
and  recharge  C  with  the  opposite  polarity.  Note,  if  we  drop  the  factor 
(l-f)"^  in  (16),  then  W.^p^/(%  CVp  is  the  fraction  of  stored  energy  trans¬ 
ferred  to  a  single  implosion. 


C.  Dimensionless  Equations. 

It  is  useful  in  deriving  scaling  laws  to  write  (9)  and  (10)  in 
the  dimensionless  form 


(17) 


where  the  dimensionless  radius  A,  current  i ,  and  time  x  are  defined  by 


j  ,  T  ~  ti)t  , 
0 


(19) 
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and 
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Solutions  are  insensitive  to  the  initial  thin  flash-over  layer  which 

appears  in  both  the  initial  value  of  A  and  also  in  the  initial  resistance 
through  the  small  quantity  e  =  6^1^5^(477003/0^)'^  in  Equation  (18).  We  note 
that  impact  with  the  projectile  of  radius  ap  occurs  when  A  =  3p/\>  st 
which  point  1  and  the  implosion  energy  is 


W.  1  =  % 
impl  ® 


m  a^  03^ 
g  w 


(^)' 
^dx'' . 
imp 


(22) 


and  the  efficiency  (16)  can  be  written 

L 


(dA/dx) 


'impl 


(di/dx)^ 

0 


2 

imp 


(1  +r^) 
91 _ ±- 


2  (1-f)  iln  (/) 
P 


(23) 


Numerical  Examples. 

From  (17)  -  (21)  we  see  that  A  and  i  are  functions  of  the  dimension¬ 


less  parameters 


For  given  values  of  these,  the  initial  value  (di/dx)^  can  be  found  numerically 

for  which  i  =  0  when  impact  occurs.  For  this  particular  case,  which  we 

choose  for  stability  reasons  (see  Section  IV),  there  is  no  inductively  stored 

energy  at  impact  and  the  capacitor  becomes  partially  recharged  with  the 

opposite  polarity.  It  is  essentially  the  case  for  which  the  implosion  time 

t.  ,  is  of  order  the  /^L  C  time  in  most  situations.  Once  we  have  numerical 
imp  I  0 

solutions  to  (17)  and  (18)  for  choices  of  the  above  parameters ,  we  can 
invert  (21)  and  use 


(24) 


(dr)  ,  (25) 

imp 

/I^  =  ^  „  (26) 

°  ^irnpl 

to  obtain  the  gas  density,  implosion  velocity,  etc. ,  for  various  systems. 

In  Figure  4,  we  have  plotted  the  dimensionless  impact  velocity  (<^A/dT)^.jj^p 
implosion  time  and  initial  gradient  (di/dt)^  (for  which  i^j^p  =  0), 

as  functions  of  the  dimensionless  resistance  r  for  the  two  examples 


^  =  0.5,  1.0; 

g 


10  . 


Using  (23) ,  the  fraction  of  energy  (T-'f)ein,pl  is  transferred  to  implo¬ 
sion  kinetic  energy  in  a  single  implosion  is  also  plotted  (f  is  the  fraction 
restored  to  the  capacitor  bank  and  potentially  recoverable  at  t^i^p-]). 

Figure  5  shows  the  variation  of  dimensionless  implosion  radius  A  and  current 
i  for  one  of  these  cases  in  the  zero  resistance  limit  r  =0  (this  would 
correspond  to  s-jj^pi  =  !)• 

We  See  from  Figure  4  that  the  transfer  efficiency  can  be  very  high,  up 
to  about  57%  for  the  case  in  which  the  external  circuit  inductance  Lg 
equals  the  maximum  module  inductance,  Ljj^,  and  the  resistance  r  is  negligible. 
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As  a  specific  example  consider  the  class  of  systems  for  which  Lp^/Lg  =  1, 
dw/ap  =  10,  r»0.  From  the  numerical  results  in  Figured,  and  noting  that 
the  initially  stored  electrical  energy  is 


W  =  Js  cv! 


0  1 


it  then  follows  from  (24)  -  (26)  that 


V, 


i.no*  W 


cm/sec 


<''kv 


gm/cm- 


,-5 


‘impl  =  ’»  (“kj)  W 


where  the  initial  voltage  and  stored  energy  W  are  in  kilovolts  and  kilo¬ 
joules  as  indicated  by  the  subscripts.  For  example,  if  a^  =  15  cm,  1=2  cm, 
Vq  =  2  kV,  W  =  5  kJ,  then  V.^p.,«10®  cm/sec,  =  3.6  lO"^  gm/cmS  and 
timpi  =16  ysec.  The  general  scaling  is  also  clear  from  the  above  formulas. 
In  order  to  achieve  higher  velocity  implosions  we  must  use  higher  driving 
voltages  or  larger  radii  modules  (although  an  upper  limit  to  the  aspect  ratio 
aw/ap  will  be  set  by  the  possible  growth  of  instabilities  in  the  imploding 
shell).  Lower  gas  densities  and  shorter  implosion  times  are  then  involved 
and  a  fuller  survey  of  this  scaling  will  be  given  later. 
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E.  Plasma  Temperature 

In  general  an  energy  equation  should  be  added  to  (9)  -  (12)  to  deter¬ 
mine  the  electron  temperature  required  in  the  conductivity,  a.  However, 
for  most  systems  the  largest  part  of  the  impedance  is  inductive,  so  that 
Eqs.  (9)  -  (12)  are  insensitive  to  T,  although  T  does  play  an  important  role 
in  the  annular  shell  thickness  as  discussed  in  Section  F. 


An  approximate  energy  equation  for  T  can  be  obtained  by  first  noting 
that  the  necessary  condition  for  a  coronal  equilibrium  plasma  state  is 
typically  satisfied,  i.e.,  n^  10^^,  since  and  the 

annular  gas  density,  which  is  compressed  over  the  ambient  value,  exceeds 
10^^/cm^  Further  the  implosion  time-integrated  Joule  dissipation  exceeds 
the  thermal  energy  content  of  the  gas,  so  that  an  approximate  balance 
between  Joule  heating  and  radiation  prevails,  i.e.. 


I^R  » smaller  of 


2TTa£  Aa 
4iTa£  cTg  T** 


(27) 


where  P^  is  the  coronal  equilibrium  emission  and  the  second  expression  the 
black  body  emission. 

Consider  for  example  nitrogen  (which  becomes  atomic  in  the  discharge) 
as  the  driving  gas.  Breton  et  al*^  give  coronal  equilibrium  results  for  the 
total  emission  (including  line,  continuum,  and  di electronic  recombination) 
that  can  be  fitted  so  that  the  above  expressions  become. 


2Tra£  Aa  4.10"PJq^^  Z  (T^^  -  3.4),  for  4^Tgy^lO 
4Tra£  lOV  watts 

■  .  .  14  ,  ■ 


Now  typical  values  for  I^R  are'^lO®  watts,  so  is  typically  a  few  eV. 
Further,  since  increases  due  to  compression  to  a  final  value  that  can 

be '''10^>  we  see  that  black-body  emission  is  approximately  dominant  for 
nitrogen.  Thus  for  finite  t. 


Since  Z^vl,  t/>»10'^  sec,  and  the  plasma  resistivity  dominates  neutral 
collisions,  this  results  in  a  temperature  which  is  clamped  at  about  2  eV. 

F.  Estimated  Thickness  of  the  Annular  Plasma  Shell. 

In  the  above  model  the  finite  thickness  Aa  of  the  imploded  plasma 
annulus  is  treated  as  small  compared  with  a.  This  thickness  determines  the 
pressure  pulse  and  characteristic  time  in  which  momentum  is  transmitted  to 
the  projectile,  and  a  proper  treatment  of  this  structure  necessitates  use 
of  a  detailed  computer  model. 

However,  a  rough  estimate  of  the  thickness  Aa  can  be  obtained  by 
considering  the  simple  case  of  a  ID  layer  of  plasma  propelled  magnetically 
with  a  constant  acceleration  g.  Assuming  field  diffuses  into  the  swept  up 
plasma  slowly  on  a  sound  transit  time  across  Aa,  an  approximate  pressure 
balance  prevails,  i.e. , 

d2 

^  +  nT  (1  +  Z)  =  -  nm^g  ,  (28) 
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from  which  the  density  profile  follows  as 


where 


is  the  sound  speed. 

Various  models  for  B  can  be  chosen,  such  as  a  linear  diffusion  ramp 
B  =  Bq(1  -  x/s),  0  <  X  <.  6  =  (c2t/4Tra)*^.  Essentially,  the  result  is  the 
profile  shown  in  Figure  6  in  which  n  increases  from  zero  on  the  swept-up 
field- supported  side  to  a  maximum  value  in  a  skin  depth  6,  and  then  under¬ 
goes  an  exponential-atmosphere  decrease  with  a  scale  length  c|/g.  A  charac¬ 
teristic  thickness  which  includes  these  two  regions  is 


/cH\^ 

VAira/ 


(29) 


A  crude  estimate  for  the  implosions  we  are  concerned  with  then  follows  by 
substituting  g»V?^p^/2a^  and  t«t^.^p^  in  (29).  For  example,  a  low  density 
air  implosion  with  T  =  2  eV,  Z  =  1,  =  10  ^  sec,  a^  =  15  cm, 

^impl  “  cm/sec,  gives  Aaa»2  cm.  A  point  to  note  is  that  high  temperature 

implosions  should  be  avoided  so  that  the  exponential  atmosphere  does  not 
extend  inward  too  far  ahead  of  the  swept-up  gas  and  thus  into  the  path  of 
the  projectile  (see  Figure  6). 
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III.  COUPLING  OF  IMPLOSION  ENERGY  TO  PROPULSION  OF  THE  PROJECTILE 

In  Figure  7  we  show  imploding  plasma  incident  with  velocity  +  V?|^pl)^ 
in  the  projectile  rest  frame,  where  U  is  the  projectile  axial  velocity  in 
the  accelerator  frame,  the  radial  plasma  velocity,  and  9  the  projectile 

cone  half  angle  along  its  impact  section.  After  impact  with  the  surface  the 
plasma  undergoes  reflection  (with  reduced  momentum  and  energy)  at  an  angle 
between  the  two  extremes  of  specular  reflection  (path  1)  or  normal  reflection 
(path  2).  Since  the  exact  reflection  vector  depends  on  details  of  the 
coupling  processes  at  the  projectile  (such  as  viscosity,  radiation,  and 
conduction  in  the  boundary  layer  between  the  solid  surface  and  the  outgoing 
shock),  we  will  confine  our  discussion  to  the  above  two  limiting  cases. 

Consider  a  series  of  modules  each  of  which  impart  a  small  velocity 
increment  dU  to  a  projectile  of  mass  M  by  imploding  a  gas  mass  dm  with 
velocity  V^.^pl .  The  projectile  energy  gain  and  implosion  energy  are 


dW„^„.  =  MU  dU 
proj 


"“impl  ' 


(30) 


Now  suppose  the  incoming  plasma  momentum  is  absorbed  and  that  a  mass  of 
gas  (dm  +  dM)  is  "re-emitted"  with  a  reduced  velocity 
an  average  angle  <t>  to  the  surface  normal  as  shown  in  Figure  7.  Momentum 
conservation  then  gives 


MdU  =  dm  0  ♦  '^(0*  +  [1  +  (5)  ^Sln  (e  +  ♦)}  ,  (31) 
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where  we  have  allowed  for  loss  of  a  fraction  (1  -  a)  of  the  implosion  kinetic 
energy  to  sinks  such  as  radiation,  ablation,  and  internal  motion  of  the 
projectile,  and  the  extra  mass  dM  re-emitted  is  due  to  ablation,  t.e. ,  the 
plasma  re-emission  velocity  was  determined  by 

a  %  dm  (U^  +  V?^p^)  =  %  (dm  +  dM) 

Note,  (dM/dm)^l^^  <  1  corresponds  to  a  mainly  implosion  momentum  drive,  and 
(dM/dm)gl^^>  1  to  a  mainly  rocket  reaction  drive  in  which  the  implosion  pro¬ 
vides  energy  but  not  much  momentum.  The  average  "re-emission"  angle  (j)  lies 
between  two  limiting  extremes: 

=  +  specul  ar  case  1 , 


(j)  =  0,  normal  emission  case  2. 


The  propulsion  energy  efficiency  can  now  be  written 


"  ^^-mpl 


2  {-  x^  +  gx  (1  +  x^)'®}  , 


where  for  the  two  limiting  cases 


3  =  VcT  1  +  (©  cos  (\p  -  29),  case  1, 

L  ^  ablJ 


3  =  vC  1  +  (^)  sine,  case  2 


X  =  Y  »  tan  lij  =  1 
^impl 


The  inequality  X  <  cot  0,  i.e. , 
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>  U  tan  6  , 


(35) 


V 


impl 


must  also  be  satisfied  in  order  for  the  gas  to  be  able  to  impact  on  the 
projectile  surface  before  it  receeds. 


Now  we  can  set  in  case  1  of  Equation  33,  since  ^ 

for  most  systems  of  interest.  With  this  replacement  the  plot  of  as 

a  function  of  x  has  two  forms,  as  shown  schematically  in  Figure  8,  depending 
on  whether  8  is  <  or  >1. 


For  cases  in  which  8<1*  the  efficiency  epj,Qp  has  a  maximum  value 
which  follows  from  (32)  as 

"  ''  ■  ’  (36) 


and  occurs  at 


(37) 


The  principal  matching  requirement  between  modules  and  projectile  needed 
to  maintain  efficient  acceleration  along  the  device  is  thus  to  match  the 
implosion  velocity  so  that  V^.^^p-j  =  U/Xn^,  provided  also  that  (35)  is  satis¬ 
fied.  For  example,  for  an  implosion  drive  situation  (or  mixed  drive)  in 
which  half  the  incident  kinetic  energy  is  lost  i.e.,  a  =0.5,  and  for 
e  =  tt/4,  we  find  e  (max) » (0.25  to  0.12)  [l  +  (dM/dm)]  ,  and  occurs  at 

P  I  W 

V.  „,»a  few  times  U. 
irnpl 

Note  also  that  in  the  ablation  drive  ^ Vijjjpi»U,  (dM/dm)gj^^»  1 , 
case  2  of  Equation  33,  we  recover  the  familiar  rocket  equation  result 

^  "  Wss  ® 
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IV,  PROJECTILE  STABILITY.  SHAPE,  AND  IMPLOSION  REQUIREMENTS 

Various  tapered  projectile  shapes  (see  Figure  3)  can  be  chosen  subject 
to  the  requirements  that  a  forward  force  results  when  the  discharge  implodes 
against  its  surface,  and  that  a  stable  progression  of  the  projectile  occurs 
as  it  advances  from  module  to  module  along  the  accelerator.  The  projectile 
must  be  stable  both  against  tumbling  and  against  cumulative  radial  drift 
perpendicular  to  the  accelerator  axis. 

To  avoid  tumbling  it  is  important  for  the  center  of  mass  of  the  projec¬ 
tile  to  be  located  back  behind  the  center  of  thrust  of  the  propelling  force. 
This  can  be  accomplished  by,  for  example,  hollowing  the  projectile  nose  or 
adding  a  slightly  extended  tail  section  as  shown  in  Figure  3.  Such  a  tail 
section  could  also  serve  to  provide  the  trailing  section  of  the  discharge 
with  a  hard  core  so  that  kink  instability  cannot  occur  as  the  main  tapered 
section  advances  into  the  next  electrode  gap  (although  since  the  magnetic 
fields  are  weak  at  impact,  kinking  is  not  expected  to  be  a  problem).  Small 
thick  fins  or  grooves  in  the  trailing  section  could  also  be  used  for  spin. 

For  accelerators  consisting  of  a  large  number  of  modules  it  will  also 
be  necessary  to  ensure  that  the  projectile  does  not  deviate  appreciably 
from  the  accelerator  axis  in  the  radial  direction.  In  Figure  9  we  show 
an  annular  implosion  converging  on  to  a  cylindrical  projectile  which  has 
become  displaced  a  distance  dr^  from  the  implosion  axis.  The  total  perpen¬ 
dicular  momentum  incident  on  a  unit  length  of  the  projectile  is  thus 

dPj.  =  1;;:  i  de  {dr_L  (1  +  sine)  ^  -  P(a)}  sin9  , 

where  P/2it  is  the  radially  directed  momentum  per  radian  of  the  plasma 
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annulus  (as  in  Eq.  (5)).  If  we  assume  that  the  re-emitted  radial  momentum 
per  radian  is  an  increasing  function  of  the  incident  momentum,  the  require¬ 
ment  for  stability  against  small  displacements  is  that 


as  shown  in  Figure  10  (note  P  is  negative).  Under  these  circumstances  the 
side  of  the  projectile  farthest  from  the  axis  experiences  a  stronger  impact 
than  the  side  closer  to  the  axis,  so  that  the  projectile  experiences  a 
radial  "impulse  well"  as  it  advances  along  the  accelerator.  The  radial  shape 
of  this  impulse  well  for  finite  amplitude  displacements,  and  assuming  various 
models  for  P( a),  has  been  analyzed  by  T.  G.  Northrop. 


In  going  from  equation  (4)  to  (9)  we  assumed  that  the  annular  shell 


thickness,  Aa,  was  sufficiently  small  that  the  thermal  pressure  term  could 
be  neglected.  However,  this  term  becomes  important  near  the  end  of  the 
implosion  so  that  (9)  should  more  precisely  be  written 


(39) 


Thus,  if  we  arrange  for  the  current  I  to  become  approximately  zero  at  impact 
(which  was  the  numerical  example  discussed  in  Section  2D),  we  see  that  the 
thermal  pressure  term  remains  in  (39)  and  causes  a  deceleration  just  prior 
to  impact,  thereby  satisfying  the  stability  requirement  (38)  (see  Figure  10  ). 


21 


Finally,  we  note  that  other  approaches  to  stability  are  also  possible 
for  long  multiple  module  devices.  For  example,  radially  offset  corrector 
modules  could  be  spaced  at  intervals  along  the  accelerator  and  provide  the 
required  re-centering  impulses,  as  computed  by  a  logical  circuit  from  sensor 
input  determining  the  trajectory  of  the  projectile,  in  which  case  matched 
implosions  might  not  be  necessary. 

V-  CONCLUDING  COMMENTS  AND  APPLICATIONS 

Development  of  the  accelerator  technology  described  here  appears  to  have 
the  advantage  that  in  its  initial  stage  it  could  be  experimentally  explored 
on  a  fairly  small  scale  (first  using  a  single  module  and  kilovolt  capacitors 
with  a  few  kilojoules).  Parameters  influencing  implosion  stability  and 
reproducibility,  timing,  the  dynamic  material  properties  of  projectiles, 
ablation,  possible  multiple  module  implosion  merging  and  smoothing,  etc., 
are  examples  of  the  areas  expected  to  play  an  important  role  in  attempting  to 
develop  long  length  accelerators. 

A  fairly  wide  range  of  accelerator  sizes  are  conceivable,  depending  on 
the  applications  intended.  These  range  from  basic  impact  studies  or  switching, 
to  more  advanced  concepts  involving  mass  accelerators.  For  example,  large 
projectiles  (kilograms)  could  be  accelerated  up  to  velocities  of  about  20 
km/sec  and  launched  as  space  probes  through  the  atmosphere  from  a  high 
altitude  site  directly  into  earth  orbit  or  into  space  (also  perhaps  for 
radioactive  waste  removal).  For  such  purposes  a  projectile  which  is  both 
compatible  with  the  accelerator  (CM  behind  the  center  of  thrust)  and  with 
aerodynamic  requirements  (CM  ahead  of  center  of  drag)  must  be  devised,  such 
as  the  two-component  projectile  shown  in  Figure  11a. 

A  second  potential  application,  this  time  involving  projectiles  of 
relatively  small  mass,  is  impact  fusion.^  This,  however,  requires  very  long 


22 


accelerators  (several  100  meters)  since  high  velocities  in  the  100  km/sec 
range  may  be  required.  It  is  interesting  to  note  that  for  ignition  of  a 
spherical  thermonuclear  pellet  the  outer  shell  must  be  driven  inward  (for 
example  with  a  laser  pulse)  with  an  implosion  velocity  of  about  200  km/sec. 
The  radius  over  which  these  implosion  velocities  can  be  reached  is  limited 
by  Rayleigh  Taylor  instability  for  large  aspect  ratio  hollow  pellets.  One 
could  view  a  conical  plug  of  fusion  fuel  embedded  in  the  nose  of  a  projec¬ 
tile  accelerated  in  a  mass  accelerator  as  a  tamped  segment  of  such  a  thermo¬ 
nuclear  pellet,  in  which  the  fuel  is  driven  inward  with  the  required  velocity 
upon  impact  with  a  target  plate  (Figure  lib).  Rayleigh  Taylor  instability 
is  avoided  as  the  "outer  shell"  velocity  relative  to  the  fuel  is  built  up 
in  the  accelerator,  because  the  projectile  and  target  plate  remain  solid 
prior  to  impact.  Simple  models  indicate  that  large  projectile  kinetic 
energies  (>  10  MJ)  would  be  required  for  ignition  in  this  particular  arrange¬ 
ment  due  to  the  limited  compressions  attainable  in  such  conical  implosions. 
However,  the  potential  advantages  of  impact  fusion  are  that  large  total 
ignition  energies  appear  to  be  accessible  with  high  accelerator  efficiency, 
and  that  the  final  transport,  aiming,  focussing,  and  small  entrance  hole 
requirements  in  the  reactor  vassal  (which  for  example  present  difficulties 
for  particle  beam  ignitors)  are  eliminated. 

Note  also  that  considerable  gains  occur  if  ignition  could  be  achieved  at 
lower  impact  velocities,  through  for  example  the  use  of  more  complex  velocity 
multiplying  projectile  or  target  structures  such  as  those  investigated  by 
Kal  iski 's^  group  (^3.10^  DD  neutrons  have  been  produced  by  purely  chemical 
explosive  means  in  these  experiments) .  The  accelerator  length  scales  as  U^, 
so  that  for  a  given  projectile  kinetic  energy  a  heavier  projectile  (>10  gm) 
could  be  used  at  lower  velocities,  which  is  also  needed  in  order  to  survive 
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the  estimated '«0.1  cm  of  surface  projectile  ablation  that  would  occur  in 
an  accelerator  in  the  100  m  length  range. 
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Annular  discharge  module  with  a  driving  voltage  in  the  1  to  tens^of  kilovolts 
range.  Power  compression  occurs  when  the  plasma  implosion  time  is  longer  than 
the  projectile  transit  time  through  the  module.  This  particular  module  is 
shown  with  a  slight  insulator  tilt  to  produce  tapered  implosions.  The 
simpler  straight  insulator  system  is  shown  in  Figure  2. 


Mass  accelerator  using  imploded  discharges  (MAID).  Assuming  the  gas  mass  m^  in  the  accelerator  is  sequentially 
swept  up  and  imploded  with  velocity  against  a  projectile  of  cone  half  angle  8  and  mass  M,  the  simple 

specular  reflection  limit  gives  U^^mg  V^^pl  sin  26 /M  for  the  final  projectile  velocity.  For  more  general 
cases  including  the  ablated  mass  reaction  see  Section  III.  More  closely  packed  multiple  module  systems  can 


Figure  3- 


Projectiles  of  either  metal  or  insulator  with  nose  cavities  which  serve  to 
move  the  center  of  mass  back  behind  the  center  of  thrust  in  the  accelerator 
for  stability.  Small  thick  oblique  tail  fins  or  grooves  could  also  be  used 
to  spin  the  projectile  if  needed.  Inner  structures  can  also  be  designed  to 
minimize  shock  convergence  onto  the  projectile  axis. 
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Plots  of  the  dimensionless  impact  velocity,  (dA/dx)  -^,  plasma  implosion  time  ,  imtial 

current  gradient  (di/dT)^,  and  implosion  efficiency,  for  two  matched  cases  (i.e. ,  cases  for y 
which  the  current  at  impact  is  zero)  as  functions  of  r  =  2S./ca Higher  efficiencies 
than  the  above  simple  circuit  model  cases  (Eqs.  10,18)  may  be  achievable  by  use  of  pulse 
forming  transmission  lines  to  reduce  the  wave  reflected  from  the  module  back  into  the  capacitc 


Plot  of  the  dimensionless  radius  A  =  a/a^  and  current  1=1/1^  as  a 
function  of  x  =  ut  for  an  implosion  in  which  electrical  resistance  of 
the  discharge  is  neglected. 


Figure  .6 

Schematic  plot  of  the  variation  of  plasma  density  n  and  magnetic 
field  for  a  layer  of  plasma  undergoing  an  acceleration  g. 


Variation  of  the  propulsion  efficiency,-  (Equation  32),  with 

‘J/^impr 


Impact  of  imploding  plasma  annulus  with  a  projectile  that  has 
become  displaced  a  distance  dr  from  the  accelerator  axis. 


Figure  10 


A  matched  implosion,  i(T^.|^p)  =  0,  showing  the  final  decrease  in  the 
implosion  momentum  due  to  plasma  thermal  pressure  forces  as  required 
to  provide  a  stabilizing  restoring  impulse  for  projectile  radial 
displacements.  Although  the  above  radially  integrated  momentum 
decreases,  the  inner  edge  of  the  annulus  would  continue  to  accelerate 
inward  and  the  outer  radial  region  to  decelerate  near  the  end  of  the 
implosion. 


Two  examples  of  complex  projectile  structures.  Case  (a)  is  an 
accelerated  structure  in  which  the  outer  jacket  is  stripped  off 
by  air  resistance  when  it  exits  the  accelerator  and  the  inner 
aerodynamical ly  suitable  projectile  proceeds  on.  Case  (b)  is 
an  i 1 1 ustration  of  a  projecti 1 e  whi ch  gi ves  ri se  to  a  tamped 
conical  compression  of  DT  thermonuclear  fuel  on  impact  with  a 
target  plate.  More  complex  velocity-multiplying  structures^ 
would  also  be  possible.  _____ 


